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We investigate the observational viability of a class of cosmological models in which the vacuum 
energy density decays linearly with the Hubble parameter, resulting in a production of cold dark 
matter particles at late times. Similarly to the flat ACDM case, there is only one free parameter to 
be adjusted by the data in this class of A(t)CDM scenarios, namely, the matter density parameter. 
To perform our analysis we use three of the most recent SNe la compilation sets (Union2, SDSS and 
Constitution) along with the current measurements of distance to the BAO peaks at z = 0.2 and 
2 = 0.35 and the position of the first acoustic peak of the CMB power spectrum. We show that in 
terms of x 2 statistics both models provide good fits to the data and similar results. A quantitative 
analysis discussing the differences in parameter estimation due to SNe light-curve fitting methods 
(SALT2 and MLCS2k2) is studied using the current SDSS and Constitution SNe la compilations. 
A matter power spectrum analysis using the 2dFGRS is also performed, providing a very good 
concordance with the constraints from the SDSS and Constitution MLCS2k2 data. 

PACS numbers: 98.80.Es, 95.35.+d, 98.62. Sb 
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I. INTRODUCTION 

The discovery of the late-time cosmic acceleration 
poses one of the greatest challenges theoretical physics 
has ever faced. In principle, this phenomenon may be 
the result of unknown physical processes involving either 
modifications of gravitation theory or the existence of 
new fields in high energy physics. In the context of Ein- 
stein's general theory of relativity this means that some 
sort of dark energy, constant or that varies slowly with 
time and space, dominates the current composition of the 
cosmos (see, e.g., [l| for recent reviews). 

Among the many possibilities to describe this exotic 
component, the simplest way is possibly by means of a 
cosmological constant A, which acts on the Einstein field 
equations as an isotropic and homogeneous source with 
p\ = — pa- Historically, A has been identified with the 
energy density of the quantum vacuum, but any attempt 
to obtain such contribution in quantum field theories in 
curved space-time leads to results many orders of magni- 
tude above the observed value (|/>a | — 10~ 10 erg/cm 3 ), 
and even when an appropriate renormalization procedure 
is used, the resulting value usually disagrees with the ac- 
tual one [2|]. This is the case, for example, of the vacuum 
density of free conformal fields in de Sitter space-time, 
which is shown to be of the order of H 4 , where H is 
the Hubble parameter [3J]. Although such a result can 
be used to construct non-singular scenarios in the high- 
energy limit [J], at late-times it leads to a very tiny A. 
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A possible way out of this problem is to consider the 
contribution of interacting fields in the low-energy limit, 
which is a more realistic situation. In this sense, some au- 
thors have argued that the QCD chiral phase transition 
originates a vacuum energy density scaling as A « m 3 H, 
where m « 150 MeV is the energy scale of the QCD vac- 
uum transition [5j. By taking A ~ H 2 in the present 
Universe, we obtain A ~ to 6 , which coincides with the 
current order of magnitude of the cosmological term. 
Clearly, no theoretical result in this field is free of some 
assumptions, and a number of other proposals for a pos- 
sible time dependence of A can be found in the literature 
(see, e.g., (g). In particular, in the case of a vacuum den- 
sity decaying linearly with H , one can show that the cor- 
responding cosmological model is indistinguishable from 
the ACDM model during the phase dominated by radi- 
ation, but it departs from the standard scenario during 
the matter epoch, in part owing to the production of dark 
matter associated to the vacuum decay [8| . 

In a previous communication [9j, we have tested the 
observational viability of this class of A(i)CDM scenar- 
ios discussed above from a joint analysis involving type 
la supernovae data (Supernova Legacy Survey Collabo- 
ration [lfj) along with the position of the first acous- 
tic peak in the CMB spectrum of anisotropics (Zi) from 
WMAP [HI, and the SDSS observed distance to the bary- 
onic acoustic oscillation (BAO) peak at z = 0.35 [12 ] . 
The goal of the present contribution is to update the re- 
sults of Ref. Q by using three of the most recent SNe 
la data sets, namely, the updated compilation of Union 
sample (Union2) [l|], the nearby + SDSS + ESSENCE 
+ SNLS + Hubble Space Telescope (HST) set of 288 SNe 
la discussed in Ref. [1J] (throughout this paper we refer 
to this set as SDSS compilation) and the Constitution 
set of 397 SNe la [28( ■ We consider two sub-samples of 
these two latter compilations that use SALT2 [15( and 
MLCS2k2 [3 SNe la light-curve fitting methods. Along 



with the SNe la data, and to improve the bounds on the 
free parameter of the model we use, besides the current 
value of h and the BAO peak at z = 0.35, the distance 
to the BAO peak at z = 0.2, as recently discussed in 
Ref. [17[. For the sake of comparison we also perform 
the same analyses for the standard ACDM model. 

In order to verify the consistence between background 
and pcrturbative tests, we also perform a joint analysis 
including the observed matter power spectrum from the 
2dFGRS. In the case of SDSS (MLCS2k2) and Constitu- 
tion (MLCS2k2 17) data, we obtain a very good concor- 
dance with A(t) model. Since MLCS2k2 depends weakly 
on standard cosmology as compared to SALT2, it is ac- 
tually the most appropriate to test non-standard models, 
which suggests that the obtained concordance is robust. 

This paper is organized as follows. In Sec. II we dis- 
cuss the basic equations and some features of the vacuum 
decay model with A ex H . Observational tests involving 
SNe la, BAO, CMB and LSS data are performed in Sec. 
III. In Sec. IV we present our main results. We also 
discuss quantitatively the influence of two SNe la light- 
curve fitting methods on parameter estimates using the 
current SDSS and Constitution compilation sets. We end 
by summarizing our main conclusion in Sec. V. 



II. BASICS OF THE A(i) MODEL 



which coincides with the Einstcin-dc Sitter solution in the 
limit of early times and at late times tends asymptotically 
to the de Sitter space-time. 

The corresponding Hubble parameter as a function of 
the redshift is given by @, [l8j 



H(z) 



H a 



i-n m + a m {i + z) 



3/2 



(5) 



where fl m and Ho = lOOh km/s/Mpc are, respectively, 
the current values of the matter density and Hubble pa- 
rameters. Note that, similarly to the flat ACDM sce- 
nario, these are the only two parameters of the model 
to be adjusted by the data. Note also that this class of 
A(t)CDM cosmologies does not contain the ACDM sce- 
nario as a particular case, which means that in principle 
they may be obscrvationally distinguishable. 



III. OBSERVATIONAL TESTS 

In this Section we investigate some observational con- 
straints on both models from five different SNe la data 
sets along with the current value of l\ from WMAP, the 
distance to the BAO peak at z = 0.2 and z = 0.35 and 
data of Large Scale Structure (LSS) distribution from 
2dFGRS. 



In units where i&irG) 1 / 2 = c = 1, the Friedmann equa- 
tions for a spatially flat universe can be written as 



p T = 3H 2 



p T + 3H(p T +p T ) = , 



(1) 



(2) 



where pr and pr are the total energy density and pres- 
sure, respectively, and the second equation expresses the 
conservation of the total energy. By taking the cosmic 
fluid as composed by dust matter of density p m and a 
time-dependent cosmological term with p\ = —A, we 
have 



l>n 



iHp n 



-A 



(3) 



which means that the decay of the vacuum energy den- 
sity is concomitant with matter production. In order to 
avoid problems with primordial nucleosynthesis and the 
spectrum of CMB, we must postulate that no baryonic 
matter or radiation is produced in this process. There- 
fore, the present model may be considered as a particular 
case of models with interaction only in the dark sector 
(see, e.g., [7| and Rcfs. therein for other examples of 
scenarios with interaction in the dark sector). 

With the ansatz A = aH, where a is a constant of the 
order of m 3 , the solution of the above equations is given 

by a 



A. SNe la data 

We use in our analyses five of the most recent SNe 
la compilations available, namely, the Union2 sample of 
Ref. [lj], the two com pila tions of the SDSS collabora- 
tions discussed in Ref. [lj] and two of the Constitution 
compilations considered in [28I] 1 . 

The Union2 sample 2 is an update of the original Union 
compilation. It comprises 557 data points including re- 
cent large samples from other surveys and uses SALT2 
for SN la light-curve fitting. The SDSS compilation 3 of 
288 SNe la uses both SALT2 and MLCS2k2 light-curve 
fitters 4 and is distributed in redshift interval 0.02 < z < 
1.55. This sample, along with the so-called CMB/BAO 
ratio, has been used in Ref. [191 ] to perform a comparative 
analysis involving a number of nonstandard cosmologi- 
cal scenarios (we refer the reader to this latter reference 



a(t) =C[exp(cr£/2)-l] 1 



(4) 



1 In |28j | the name Constitution refers to data calibrated with 
SALT. Here we will use it to refer to any compilation presented 
there. 

2 http://www.supernova.lbl.gov/ 

3 http://www.sdss.org/ 

4 It is worth emphasizing that SALT2 fitter does not provide a 
cosmology-independent distance estimate for each SNe, since 
some parameters in the calibration process are determined in a 
simultaneous fit with cosmological parameters to the Hubble di- 
agram (see 12011 for a discussion). See also I2l]l for a discussion on 
possible systematic biases from the multicolor light curve shape 
method. 



for more on the SDSS compilation). The Constitution 
compilations comprises 397 SNe la (CfA3 + Union) cal- 
ibrated with SALT2 and two versions of the MLCS2k2 
fitter: MLCS31 and MLCS17. As discussed in [H , the 
former overestimates host-galaxy extinction, and here we 
will consider only the later version. 



B. BAO 

The second observable we consider is the distance to 
a given observed baryonic acoustic oscillation, which is 
defined as 



D v (z) 



D M {zf 



H{z) 



1/3 



(6) 



where Dm is the comoving angular-diameter distance. 
We will consider two different measurements of BAO, at 
redshifts zbao = 0.20 and zbao = 0.35, i.e., 

D v (z = 0.2) = 748 ± 57 Mpc, 



D v {z = 0.35) = 1300 ± 88 Mpc, 

where the value of Dy at z = 0.2 is obtained by combin- 
ing the results of Ref. [17| with the value of Dy(z = 0.35) 
obtained in Ref. [22j |. 

An important aspect worth emphasizing at this point 
is that we cannot use in our analysis the parameter A, 
defined as A = Dy(z)\/O m ffg/zBAO 01' since m the 
present model the production of matter leads to a dif- 
ferent dependence of the horizon at the time of matter- 
radiation equality on the present matter density [9j . In all 
the cases, the observable quantity is the ratio r s {zd) / 'Dy , 
where r s (zd) is the sound horizon at the time of drag 
epoch. Since r s (zd) is very weakly model-dependent (pro- 
vided the radiation epoch is the same as in the ACDM 
model), the distance Dy can in fact be used in our anal- 
ysis. As well known, it is also obtained from direct obser- 
vations by using the ACDM as a fiducial model [12[ , but 
such a dependence is weak, and actually this particular 
observable does not impose strong limits on the model 
parameter [j|. 



C. CMB 

The most restrictive of the background tests we will 
perform in our analysis is the position of the first peak 
of the CMB spectrum of anisotropics. The simplest way 
to use it is certainly by considering the shift parameter, 
defined as 7Z = ^/TE^Dm(z*), where z„ is the redshift of 
the last-scattering surface. But again, this is not allowed 
in our analysis [9J, since the angular scale I a is not the 



same as in the ACDM case 5 . We, therefore, use dir ectly 
the observed position of the peak, l\ = 220.8 ± 0.7 
which is related to the angular scale by [24] 
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h = l A (l-6 1 ), where <5 X = 0.267 



r 

03 



(7) 



with r = p r {zis)/p m (zi s ). Here, p r {zi s ) is the radiation 
density at the time of last scattering 6 . Properly speaking, 
these relations were obtained for the ACDM case, but 
they depend only on pre-recombination physics and can 
be applied to the present model as a good approximation. 
Nevertheless, owing to the process of matter production, 
the matter density at the time of last scattering is p m = 
iH^ 2 ll zf sl i.e., smaller than the ACDM value by a factor 
Q m , what can be seen directly from the high-z limit of 
the Hubble function ([5]). Therefore, in the present case 
we have [9| 



VL 2 



-Zls 



(8) 



leading to a relation between I a and l± different from the 
ACDM case. 

The theoretical expression for the acoustic scale is de- 
fined as the ratio 



Ia 



dz 

H[z)' 



dz 



c H{ z y 



9 n b 
4 n~. 



(9) 



-1/2 



with the sound velocity given by c s = c I 3 • 

where Qb and f2 7 stand for the present density param- 
eters of baryons and photons, respectively. The Hubble 
function used in ((9|) must take into account the radiation 
contribution fl r and is rewritten as |9| 



H(z) 
Ha 



i - n m + n m (i + z) 3/2 ] +n r (i + z) 



1/2 



(10) 



D. The matter power spectrum 

Except for scales near the horizon, vacuum perturba- 
tions are negligible (see [25|, [29[). Therefore, we sim- 
plify our analysis by considering an unperturbed time- 
dependent cosmological term. The model prediction 
(blue line) together with the 2dFGRS data points [3fj| 
are shown in the left panel of Figure 6. For compari- 
son, the concordance ACDM fitting as provided by the 
BBKS transfer function [31j is also shown (red line), as 



5 For the same reason, we cannot use in the joint analysis the ratio 
Ify/Dy, used in other analyses. 

6 As shown in Q , in the present model the redshift of last scatter- 
ing is the same as in the ACDM case. 



TABLE I: Limits to fi m (A(t)CDM). 



Test 



a, 



2 2 

Xmin X.min 



I" 



Union2(SALT2) 0.339±g;gtl 

Union2(SALT2) i ' 0.420±g;g?g 

SDSS (SALT2) 0.352±g;g£ 

SDSS (SALT2) 6 0.430±g ™| 

SDSS (MLCS2k2) 0.484±°;°?j? 

SDSS (MLCS2k2) 6 0.450to° w 

Constitution (SALT2) 0.355±g;g^ 

Constitution (SALT2) 6 0.450±g;g" 

Constitution (MLCS2k2) 0.402±°^4g 

Constitution (MLCS2k2) 6 0.450±g;g?l 



"Error bars stand for 2er 
*>+ CMB + BAO+LSS 



TABLE II: Limits to fi m (ACDM) 



554.9 


0.998 


634.5 


1.063 


256.9 


0.895 


320.8 


0.978 


249.7 


0.870 


276.3 


0.842 


343.3 


0.984 


390.2 


0.998 


397.2 


1.074 


435.3 


1.057 



Test 



■-'/-, 



2 2 / 

Xmin Xmin / ^ 



Union2(SALT2) 0.270±o.o38 

Union2(SALT2) 6 0.235 ± 0.011 

SDSS (SALT2) 0.273±°;°^ 

SDSS (SALT2) b 0.226 ± 0.012 

SDSS (MLCS2k2) 0.3991°;°^ 

SDSS (MLCS2k2) 6 0.2601°;°^ 

Constitution (SALT2) 0.282i°:o49 

Constitution (SALT2) 6 0.265±g;gil 

Constitution (MLCS2k2) 0.327±g;g|? 

Constitution (MLCS2k2)'' 0.270 ± 0.013 



552.7 
613.1 
256.0 
305.8 
249.1 
403.8 
341.7 
515.8 
397.4 
570.2 



0.994 
1.027 
0.892 
0.932 
0.868 
1.231 
0.979 
1.319 
1.074 
1.384 



"Error bars stand for 2a 
b + CMB + BAO+LSS 



well as the result for the ACDM case of the same sim- 
plified analysis we have performed with A(t)CDM (black 
line). In such analyses baryons are not included, which 
leads to an error of about 10% as compared to a more 
exact study [25J. In the case of ACDM, the best fit is 
given by f2 m sa 0.2, in agreement with [30(. 



IV. ANALYSIS AND RESULTS 



In order to derive the constraints on our model we max- 
imaze the likelihood function C ex exp (— x 2 /2), where 

X 2 = xInc + Xbao + Xcmb + Xlss takes int0 account all 
the data sets discussed above. In our analyses, the Hub- 
ble parameter H is considered as a nuisanse parameter 
so that we marginalize over it. The results are shown in 



Figs. 1-5, and Tables I and II summarize the main results 
of this paper for all SNe data sets and joint analyses. In 
order to compare the two classes of models discussed in 
this paper we also show the values of Xmin an d Xmin/" 
(y stands for degrees of freedom) for each analysis per- 
formed. 

In the figures, the left and central panels show the su- 
perposition of the lcr and 2<r confidence regions in the 
h — f2 m plane for the four tests we are considering, as 
well as the concordance ellipses of our joint analyses. In 
the case of SDSS (MLCS2k2), the interval for the matter 
density in A(t)CDM is about 10% higher than the fiat 
ACDM value. Both values, on the other hand, are higher 
than the usually accepted concordance value Ci m « 0.27 
fill ] - We also note that in both cases the value of H a 
is relatively small (h ~ 0.64) compared to current esti- 
mates [32j. Nevertheless, when we marginalize over h the 
matter density remains unchanged in both models. The 
right panels of the figures show the likelihood for TO af- 
ter marginalization. In all the analyses the estimate of 
the matter density parameter for the A(t)CDM scenario 
is considerably larger than the corresponding value for 
the ACDM model. In particular, the estimates of O m 
for the vacuum decay model is up to ~ 2<r off from the 
central value inferred by applying the virial theorem to 
cluster dynamics [26[ (see also [27|), but only ~ 0.5er off 
from the central value obtained by using other indepen- 
dent methods, such as the mean relative peculiar velocity 
measurements for pairs of galaxies [23[ . 

In terms of parameter estimation we note that very 
similar results are found for both Union2 and SDSS com- 
pilation with SALT2. An interesting aspect worth em- 
phasizing is realized when a comparison between the 
analyses of Figs. 2 and 3 is made. Note that, al- 
though the only difference between these two analyses 
is the light-curve fitter used in the SNe calibration, the 
estimated values of the cosmological parameter £7 m are 
considerably modified from analysis to analysis. For in- 
stance, by comparing only the SNe la bounds on fl m for 
the standard ACDM model, we observe that the esti- 
mated values of the matter density parameter increases 
- 47%, from 0.27 to 0.40 (see Table II). The same can 
be said about the Constitution sample calibrated with 
SALT2 and MLCS17. Similar conclusions arc also drawn 
for the class of vacuum decay models studied here. In 
this case, estimates of the matter density parameter in- 
crease ~ 38% for the SDSS SNe la data with MLCS2k2 
light-curve fitter, in better accordance with the indepen- 
dent estimate based on large scale structure distribution 
|25| . It is worth observing that these estimated values of 
fl m obtained from SDSS (MLCS2k2) analyses for both 
models are relatively large, being ~ 2<r off from the cen- 
tral value obtained using the cluster dynamics methods 
discussed earlier. 

In all the cases above we have considered spatially flat 
scenarios, usually predicted by inflation. For complct- 
ness, we can relax this prior and include spatial curvature 
in both models. In the right panel of Figure 6 we show 
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FIG. 1: a) Superposition of the confidence regions in the h x Q m plane from Union2 SNe la, CMB, BAO and LSS data sets 
discussed in the text, for the A(t)CDM model. The blue elipse shows the joint 2cr confidence region, b) The same as in Panel 
(a) for the ACDM model, c) The likelihood £ for Q m arising from Union2 SNe la sample for A(t)CDM (dot-dashed blue 
curve) and ACDM (dot-dashed black curve) scenarios. The solid lines correspond to the joint analysis including CMB, BAO 
and LSS data sets. 
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FIG. 2: The same as in Figure 1 for SDSS SNe la data with SALT2 light-curve fitter. 



the results for the SDSS (MLCS2k2) data set ]3$ . We 
can see that the spatially flat universe is in accordance 
with data within la confidence level in both models. 



FINAL REMARKS 



Although in good agreement with current observations, 
a residual A term exacerbates the well known cosmologi- 
cal constant problem, requiring a natural explanation for 
its small, but nonzero, value. In this paper, we have dis- 
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FIG. 3: The same as in Figure 1 for SDSS SNe la data with MLCS2k2 light-curve fitter. By comparing Figs. 2 and 3 we clearly 
see the influence of SNe la light-curve fitting on parameter estimation. 
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FIG. 4: The same as in Figure 1 for Constitution SNe la data with SALT2 light-curve fitter. 



cussed the observational viability of a class of alternative 
A(i)CDM scenarios in which the vacuum energy density 
decays linearly with the Hubble parameter H. By con- 
sidering the most recent SNe la compilations sets along 
with the current measurements of distance to the BAO 
peaks at z = 0.2 and z = 0.35 and the position of the first 
acoustic peak of the CMB power spectrum, we show that 
in terms of x 2 statistics both models provide good fits to 
the data with the very same number of free parameters 



(see Tables I and II). A recurrent aspect in all the analy- 
ses performed is the larger estimates of J7 m for A(i)CDM 
relative to the values obtained in the ACDM context. 
When compared with some fl™ estimates from indepen- 
dent observational methods [2J, [2(1 [27| , these values may 
be off from ~ 0.5er up to ~ 2<r. The same conclusion also 
arises in the ACDM context for the tests including SDSS 
(MLCS2k2). 
We have also discussed quantitatively differences in pa- 
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FIG. 5: The same as in Figure 1 for Constitution SNe la data with MLCS17 light-curve fitter. 
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FIG. 6: Left. The theoretical power spectrum for the A(t)CDM model (blue line), for the simplified (without barions) ACDM 
model (black line) and for the BBKS transfer function (red line). The data points are taken from [30] (2dFGRS). Right. The 
confidence regions in the Qa X f2 m plane for the ACDM and A(t)CDM models with spatial curvature, when we use the SDSS 
(MLCS2k2) SNe la compilation. The lines labled k — and go = correspond, respectively, to spatially flat and uniformely 
expanding universes. 



rameter estimates between SALT2 and MLCS2k2 light- 
curve fitters using the current SDSS and Constitution 
SNe la compilations. We have shown that for both mod- 
els these differences may reach up to ~ 47% for the mat- 
ter density parameter. These results reinforce the need 
and importance for a better understanding of unresolved 
systematic effects on SNe la calibration, as has been dis- 



cussed by other authors (see, e.g., (l9l|). 

Finally, from a joint analysis which includes the ob- 
served matter power spectrum, it is clear that a good 
concordance is achived for the A(t)CDM model when we 
use the SDSS (MLCS2k2) and Constitution (MLCS17) 
compilations of SNe la. Since the MLCS2k2 fitter de- 
pends weakly on standard cosmology, it is generally con- 



sidcrcd the most appropriate to test alternative models, 
which means that the concordance obtained here can be 
considered robust. Actually, an inspection of panels (c) 
of Figures 1-5 shows that, in the case of MLCS2k2 fit- 
ter, the A(i) model presents better concordance than the 
standard model. 



ful for the support of CNPq grants 304569/2007-0 and 
481784/2008-0. SC is thankful for the CNPq grant 
305133/2008-0. 



Acknowledgements 

The authors acknowledge CNPq - Brazil for the grants 
under which this work was carried out. JSA is grate- 



[1] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. 
D9, 373 (2000); P. J. E. Peebles and B. Ratra Rev. 
Mod. Phys. 75, 559 (2003); T. Padmanabhan, Phys. 
Rept. 380, 235 (2003); E. J. Copeland, M. Sami 
and S. Tsujikawa, Int. J. Mod. Phys. D15, 1753 
(2006) ; J. S. Alcaniz, B raz. J. Phys. 36, 1109 (2006) 
[arXiv: astro^p h/060863 1 . 

[2] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 

[3] L. H. Ford, Phys. Rev. Dll, 3370 (1975); J. S. Dowker 
and R. Critchley, Phys. Rev. D13, 3224 (1976); P. C. W. 
Davies, Phys. Lett. B68, 402 (1977); A. A. Starobinsky, 
Phys. Lett. B91, 99 (1980). 

[4] S. Carneiro and R. Tavakol, Gen. Rel. Grav. 41, 2287 
(2009). 

[5] R. Schutzhold, Phys. Rev. Lett. 89, 081302 (2002); F. R. 
Klinkhamer and G. E. Volovik, Phys. Rev. D 79, 063527 
(2009); F. R. Urban and A. R. Zhitnitsky, Phys. Rev. 
D80, 063001 (2009); Phys. Lett. B688, 9 (2009); Nucl. 
Phys. B835, 135 (2010); N. Ohta, arXiv: 1010.1339. See 
also S. Banerjee et al, Phys. Lett. B611, 27 (2005). 

[6] M. Ozer and M. O. Taha, Phys. Lett. B171, 363 (1986); 
Nucl. Phys. B287, 776 (1987); O. Bertolami, Nuovo Ci- 
mento Soc. Ital. Fis. B93, 36 (1986); K. Freese et al, 
Nucl. Phys. B287, 797 (1987); W. Chen and Y-S. Wu, 
Phys. Rev. D41, 695 (1990); M. S. Berman, Phys. Rev. 
43, 1075 (1991); D. Pavon, Phys. Rev. D43, 375 (1991); 
J. C. Carvalho, J. A. S. Lima and I. Waga, Phys. Rev. 
D46 2404 (1992); A. I. Arbab and A. M. M. Abdel- 
Rahman, Phys. Rev. D50, 7725 (1994); J. A. S. Lima 
and J. M. F. Maia, Phys. Rev D49, 5597 (1994); J. M. 
Overduin and F. I. Cooperstock, Phys. Rev. D58, 043506 
(1998); J. M. Overduin, Astrophys. J. 517, LI (1999); J. 
Shapiro and J. Sola, Phys. Lett B475, 236 (1999); M. V. 
John and K.B. Joseph, Phys. Rev. D61, 087304 (2000); 
O. Bertolami and P. J. Martins, Phys. Rev. D61, 064007 
(2000); L. Amendola, Phys. Rev. D62 043511 (2000); L. 
Amendola and D. Tocchini-Valentini, Phys. Rev. D 64, 
043509 (2001); R. G. Vishwakarma, Gen. Rel. Grav. 33, 
1973 (2001); A. S. Al-Rawaf, Mod. Phys. Lett. A14, 633 
(2001); M. K. Mak, J. A. Belinchon, and T. Harko, IJMP 
Dll, 1265 (2002); I. L. Shapiro and J. Sola, JHEP 0202, 
006 (2002); W. Zimdahl and D. Pavon, Gen. Rel. Grav. 
35, 413 (2003); M. R. Mbonye, IJMP A18, 811 (2003); 
J. S. Alc aniz and J. M. F. M aia, Phys. Rev. D67, 043502 
(2003). [astro-ph/0212510l ; I. L. Shapiro, J. Sola, C. 
Espana-Bonet, and P. Ruiz-Lapuente, Phys. Lett. B574, 



149 (2003); R. Horvat, Phys. Rev. D70, 087301 (2004); 

C. Espana-Bonet et al. JCAP 0402, 006 (2004); P. Wang 
and X. Meng, Class. Quant. Grav. 22, 283 (2005); S. 
Carneiro and J. A. S. Lima, IJMP A20 2465 (2005); I. 
L. Shapiro, J. Sola, and H. Stefancic, JCAP 0501, 012 
(2005); E. Elizalde, S. Nojiri, S. D. Odintsov, and P. 
Wang, Phys. Rev. D71, 103504 (2005); R. Aldrovandi, 
J. P. Beltran Almeida, and J. G. Pereira, Grav. & Cos- 
mol. 11, 277 (2005); F. Bauer, Class. Quant. Grav. 22, 
3533 (2005); J. S. Alcaniz and J. A. S. Lima, Phys. Rev. 
D 72, 063516 (2005); J. D. Barrow and T. Clifton, Phys. 
Rev. D 73, 103520 (2006); B. Wang, C. Y. Lin, and E. 
Abdalla, Phys. Lett. B637, 357 (2006); J. Sola and H. 
Stefancic, Mod. Phys. Lett. A21, 479 (2006); A. E. Mon- 
tenegro Jr. and S. Carneiro, Class. Quant. Grav. 24, 313 
(2007); F. E. M. Costa, J. S. Alcaniz and J. M. F. Maia, 
Phys. Rev. D 77, 083516 (2008); G. Caldera-Cabral, 
R. Maartens and L. A. Urena-Lopez, Phys. Rev. D 79, 
063518 (2009); 

[7] J. F. Jesus, R. C. Santos, J. S. Alcaniz and J. A. S. Lima, 
Phys. Rev. D 78, 063514 (2008); Q. Wu, Y. Gong, 
A. Wang and J. S. Alcaniz, Phys. Lett. B 659, 34 (2008); 
F. E. M. Costa, E. M. . Barboza and J. S. Alcaniz, Phys. 
Rev. D 79, 127302 (2009); F. E. M. Costa and J. S. Al- 
caniz, Phys. Rev. D 81, 043506 (2010); 

[8] H. A. Borges and S. Carneiro, Gen. Rel. Grav. 37, 1385 
(2005). 

[9] S. Carneiro, M. A. Dantas, C. Pigozzo and J. S. Alcaniz, 
Phys. Rev. D77, 083504 (2008). 
[10] P. Astier et al, Astron. Astrophys. 447, 31 (2006). 
[11] G. Hinshaw et al, Astrophys. J. Suppl. 170, 288 (2007); 

D. N. Spergel et al., Astrophys. J. Suppl. 170, 377 (2007). 
[12] D. J. Eisenstein et al, Astrophys. J. 633, 560 (2005). 
[13] M. Kowalski et al., Astrophys. J. 686, 749 (2008). 

[14] R. Kessler et al., Astrophys. J. Suppl. Ser. 185, 32 (2009). 

[15] J. Guy et al, Astron. Astrophys. 466, 11 (2007). 

[16] M. M. Phillips, Astrophys. J. 413, L105 (1993); A. G. 

Riess, W. H. Press, and R. P. Kirshner, Astrophys. J. 

438, L17 (1995); S. Jha, A. G. Riess, and R. P. Kirshner, 

Astrophys. J. 659, 122 (2007). 
[17] W. J. Percival et al, Astrophys. J. 657 51 (2007). 
[18] S. Carneiro, C. Pigozzo, H. A. Borges and J. S. Alcaniz, 

Phys. Rev. D74, 23532 (2006). 
[19] J. Sollerman et al., Astrophys. J. 703, 1374 (2009). See 

also G. R. Bengochea, arXiv: 1010.4014. 
[20] J. A. Frieman, AIP Conf. Proc. 1057, 87 (2008) 



[21] R. Kessler et al., Astrophys. J. Suppl. Ser. 185, 32 (2009) [28] 

[22] M. Tegmark et al., Phys. Rev. D74 (2006) 123507 [29] 

[23] Feldman, H. A., et al. 2003, Astrophys. J. 596, L131 
[24] W. Hu et al, Astrophys. J. 549, 669 (2001). [30] 

[25] H. A. Borges, S. Carneiro, J. C. Fabris and C. Pigozzo, [31] 

Phys. Rev. D77, 043513 (2008). 
[26] R.G. Carlberg et al., Astrophys. J. 462, 32 (1996). 
[27] A. Dekel, D. Burstein and S. D. White, "Measuring [32] 
Omega", in Critical Dialogues in Cosmology (proceed- 
ings of the Princeton 250th Anniversary conference, June [33] 
1996), ed. N. Turok (World Scientific). N. A. Bahcall and 
X. Fan, Publ. Nat. Acad. Sci. 95, 5956 (1998). 



M. Hicken et al, Astrophys. J. 700, 1097 (2009). 

W. Zimdahl, H. A. Borges, S. Carneiro, J. C. Fabris and 

W. S. Hipolito-Ricaldi, JCAP 1104, 028 (2011) . 

S. Cole et al, MNRAS 362, 505 (2005). 

J. Martin, A. Riazuelo and M. Sakellariadou, Phys. Rev. 

D61, 083518 (2000); J. M. Bardeen, J. R. Bond, N. 

Kaiser and A. S. Szalay, Astrophys. J. 304, 15 (1986). 

W. L. Freedman et al. [HST Collaboration], Astrophys. 

J. 553, 47 (2001). 

C. Pigozzo, PhD thesis, 2010. 



